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dimensions. In accordance with the prediction no qualitative dependence on the gauge 
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for becoming large in the SU(A) case. Quantitative effects are found, and analyzed in 
detail. 
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1 Introduction 

Yang-Mills theories are gauge theories of (semi-)simple Lie-algebra valued gauge fields A^, 
governed by the Lagrangian [1] 

r — _ 1 pfJ.i^a pa 

where the /"'"^ are the structure constants of the chosen gauge-algebra. The coupling 
constant g determines, in four dimensions via dimensional transmutation, in other dimen- 
sions directly, the scale of the theory. The gauge field themselves belong to the adjoint 
representation of the gauge algebra. 

It is always possible to formulate Yang-Mills theories in terms of group-valued variables. 
However, at face value all group representations for a given gauge algebra would provide 
the same dynamics on the level of the gluon fields. E. g., for the case of the su(A^) algebra, 
a representation in terms of the groups SU(A^) and S'U{N)/Zf^ must be equivalent. This is 
trivial so, as both representations differ only by group-elements, the center, which leave the 
gauge fields themselves invariant: Center transformations act as identity transformations 
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on the gluons. However, though the dynamics is necessarily the same, it is of course 
possible to construct group- valued quantities, which will not be the same in both cases, 
e. g., Polyakov loops^. Such constructions become even even more important when the 
theories are formulated on a non-trivial space-time manifold, like a torus in lattice gauge 
theory or at finite temperature in equilibrium. Anomalies due to matter fields can restrict 
the possible group representations further. In case of the standard model actually the only 
permitted group representation of the gauge algebra su(3) xsu(2) xu(l) is S(U(3) xU(2)). 
This yields the unbroken gauge group SU(3)/Z3xU(l) of the electrostrong interactions and 
the broken SU(2)/Z2 of the weak isospin [2]. 

However, the dynamics could be very different for different gauge algebras. In pertur- 
bation theory, this is evidently not the case: The different algebras only manifest themselves 
in different coefficients of the expansion, without altering qualitatively the behavior [1]. In 
particular, for any Lie-algebra Yang-Mills theories are (in four dimensions) renormalizable 
with the same number of independent renormalization constants, and are all asymptoti- 
cally free. This property is not remaining beyond perturbation theory. E. g., the order of 
the finite-temperature phase transition depends on the chosen gauge algebra, as already 
the comparison of su(2) and su(3) in terms of the groups of SU(2) and SU(3) shows [4]: 
In one case the transition is second order, in the other first order. Also the bound state 
spectrum is different due to algebraic reasons. E. g., su(2) is not having negative charge- 
conjugation-parity glueball bound states due to the existence of real representations, while 
su(3) possesses this type of bound states [5]. 

As soon as matter fields in any representation are coupled to the theory the behavior 
can even change totally. Even asymptotic freedom is no longer guaranteed if a sufficient 
number of additional matter fields in appropriate representations are added to the theory. 
This much more complicated problem will not be treated here, and only the pure Yang-Mills 
case will be studied. 

The main focus of interest here is whether the gluons show different behavior for differ- 
ing gauge algebras. This will be investigated using lattice gauge theory, by implementing 
particular gauge groups, and utilizing the gauge-dependent gluon correlation functions. 
This complements investigations of gauge-independent quantities, like the glueball spec- 
trum [6], string tensions [7], or the thermodynamic properties [8-11]. The advantage of 
investigating the gauge-dependent correlation functions is that they can be assembled to 
obtain gauge-invariant quantities, see e. g. [12-14], in analogy to perturbation theory. This 
is an approach which can be applied even were pure lattice calculations are yet restricted, 
like the chiral limit [12, 13], cold and dense matter [14-16], or cases with disparate charac- 
teristic scales [12]. Analyzing the group-dependence of the correlation functions is then an 
important input to improve the systematics of the approach and of the assumptions made. 

Since the properties of gluons are gauge-dependent this has to be done in a fixed gauge. 
This will be chosen here to be the Landau gauge. In the non-perturbative domain this is 
not a sufficiently precise definition of the gauge due to the presence of Gribov-Singer copies 

'^In case of semi-simple Lie groups, this can become even more complicated, since the algebra is just a 
product-algebra, while this is not necessarily the case for the corresponding groups [2, 3]. 
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[17, 18]. For the present purpose the minimal [19] or average-B [20] Landau gauge as the 
sub-type of non-perturbative Landau gauges will be chosen. However, since the calculation 
of large gauge-algebras becomes exceedingly expensive when it comes to gauge-fixing, only 
rather small volumes will be studied. For these volumes, a significant effect of Gribov 
copies for the gluon is not expected, and only a moderate effect for the ghost [20-22]. 

Furthermore, due to the numerical costs for the calculations, here only the case for two 
and three dimensions will be investigated, and only the propagators will be determined. In 
case of the gauge groups SU(2) and SU(3) also comparisons in four dimensions have been 
performed [23-27]. These showed no essential differences in both cases^. As this agrees 
with the findings here, it is conceivable that also the results for other gauge groups can 
be translated from lower dimensions to four dimensions, though this conjecture requires 
confirmation. The set of groups covered here is SU(2), SU(3), SU(4), SU(5), SU(6), and 

G2. 

The ladder of SU(A^) groups permits to investigate also the possibility to which ex- 
tent at fixed 't Hooft coupling X = g^N the results show a simple dependence on N with 
increasing N. Indeed, many results on gauge-invariant quantities seem to indicate such a 
behavior, see e. g. [4-6, 28-30]. All of these quantities have a genuine non-perturbative 
origin, and require confinement as well as strong interactions for the bound states to exist 
[30], and therefore indicate that the large- limit is non-perturbative. This is to be ex- 
pected, already based on Haag's theorem [31], irrespective of whether the possible range 
of possible g values, and therefore of A, is A^-dependent or not. Here, evidence will be 
provided that the gauge-dependent correlation functions are essentially independent of A^ 
at fixed A, confirming that the large- A^ limit is highly non-perturbative, but exhibiting a 
simple scaling with A^, indeed in this case a trivial one. 

The motivation to study G2 is different. Since all realizations of the g2 algebra have 
a trivial center, many arguments based on degrees of freedom sensitive to the center of 
the gauge group do not apply [11, 32]. Nonetheless, the dynamics of G2 gauge theories 
appears to be quite similar to theories with a center [10, 11, 33-37]. This observation will 
also be confirmed here. To understand this result and its implications is therefore helpful 
to identify the relevant effective degrees of freedom in gauge theories. In a sense, since g2 
has a rather different group-theoretical structure, it tests the extreme case of choosing a 
different gauge algebra. The investigation of G2 Landau-gauge propagators here extends 
previous studies [33]. 

The insensitivity on the gauge group and the non-triviality of the N 00 limit are 
actually not a surprise, and have been pointed out in functional calculations [38-40]. These 
arguments will be rehearsed in subsection 2.3. In the corresponding section 2 also other 
aspects of the continuum theory will be discussed. However, calculations using functional 
methods involve truncations, and therefore the results have to be tested. This is the main 
aim here, to test these predictions using the complementary approach of lattice gauge 
theory. The setup for these lattice calculations is given in section 3 and 4. The results, 
demonstrating indeed a qualitative gauge-group independence, are given in section 5. A 

^See, however, [27] for a differing conclusion in four dimensions. 
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short summary will be provided in section 6. 

2 The asymptotic gauge-group dependence 
2.1 Setup 

Continuum (Euclidean) Yang-Mills theory in (perturbative) Landau gauge is formulated 
as the Lagrangian 

I)^' = rt^ + ^r''^^- (2.1) 

Herein are (c'^) the (anti-)ghost field, and the index a on both fields counts the dimen- 
sionality of the adjoint representation of the gauge group. 

To complete the gauge non-perturbatively, it is necessary to chose between non- 
perturbative realizations [17, 18]. Arguments have been provided that this can be done 
by imposing conditions on certain correlation functions [20, 21, 41], e. g. on the ghost 
propagator [20, 41]. However, this issue is not yet fully settled, see [42] for a status report. 

Here, for the sake of simplicity, it will be assumed that the construction of [20, 41] is 
indeed possible. Then, a gauge choice can be made such that the ghost dressing function 
p'^Dg{p), where Dq is the ghost propagator, is infrared singular. This permits a simpler 
calculation when using functional equations below. A similar line of argumentation can be 
expected to hold also in photon-ghost gauges, where p^Dg{p) is chosen to be finite, though 
this has yet to be done explicitly. For the lattice calculations presented in section 3, it 
can be expected that the choice of non-perturbative Landau gauge is not relevant for the 
gauge-group dependence, in particular due to the limited volumes accessible. Therefore, 
the numerically much cheaper minimal or average-B Landau gauge is used [19, 20], which 
belongs to the class of photon- ghost gauges. 

One possibility to obtain from the Lagrangian (2.1) the propagators are the Dyson- 
Schwinger equations (DSEs) [12]. Keeping all color indices and the structure constants 
explicit, the Dyson-Schwinger equations are form- invariant for all gauge groups. In partic- 
ular, only the structure constants /"^'^ do appear. Additional structure constants, like the 
symmetric one 

^ahc^ being zero for most gauge groups [43], do not appear explicitly in the 
DSEs. Of course, the full Green's functions can develop such contributions. 

The explicit derivation of the DSEs is discussed in great detail elsewhere [12]. Only 
two particular limits will be of interest here. One is the far infrared limit, and the other the 
far ultraviolet limit. In the latter case, the equations generate (resummed) perturbation 
theory. 

The following discussion is presented for completeness, and to emphasize the role of 
the gauge group. It follows previous presentations of the subject [39, 40]. 
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2.2 The far ultraviolet 

There is a qualitative difference between four-dimensional and lower-dimensional systems: 
Irrespective of the gauge group neither in two, nor in three dimensions any physical renor- 
malization occurs^. 

In four dimensions the propagators behave at large momenta to leading order like 

p^D{p)'"''A^''il+iv\np)', (2.2) 

where to and 5 depend on the gauge group, but have the same sign for all gauge groups, and 
Aym is the scale generated by dimensional transmutation. In two and three dimensions, 
however, the propagators depend only polynomial on the momentum. In particular, no 
resummation occurs to leading order in g, and only starting from next-to-leading order 
such effects appear [45]. Hence, the propagator take at large momenta the form 

p^i^(p)^^'l + c^ (2.3) 
P 

in three dimensions and 

p D{p) 1 + c — 5— (2.4) 

in two dimensions, just on dimensional grounds. Ca is the adjoint Casimir of the gauge 
group, defined by 

Ca^"^ J acd jbcd 

The constants of proportionality c depend only on the underlying space-time manifold. 
In three dimensions they take the value 11/64 for the gluon propagator and 1/16 for the 
ghost propagator. In two dimensions infrared divergencies appear, which make a somewhat 
ad-hoc regularization necessary^. This already implies that non-perturbative contributions 
will appear which will provide such a regularization. It turns out that this is due to an 
infrared suppression, compared to tree-level, of the gluon propagator [47-50]. This makes 
the results convergent. 

Hence, already in perturbation theory the gluon (and ghost) propagator differ for the 
different algebras. E.g., the leading order perturbative coefficient changes by the ratio of 
the adjoint Casimir operator from gauge algebra to gauge algebra in two and three di- 
mensions. However, the momentum dependence to first order stays the same. Therefore, 
perturbatively, the gauge algebra dependence is only quantitative. Note that the pertur- 
bative expansion does not depend on the chosen group representation of the algebra. 

2.3 The far infrared 

For the sake of simplicity the argument why the gauge-algebra dependence in the non- 
perturbative domain should also be only a quantitative rather than a qualitative effect, the 



^Note that regularization is still necessary. However, a standard BPHZ-like prescription would permit 
to perform this implicitly [44]. 

* Strictly speaking, this problems occur also in three dimensions, though only at higher order [46]. I am 
grateful to David Dudal for pointing this out. 
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following discussion will be done in a scaling-type gauge [20, 41]. This section essentially 
follows the line of arguments of [40]. The situation in a photon-ghost gauge [41, 51-54], 
like the one used for the lattice calculations, will be discussed below. 

Here, only the leading infrared contributions will be retained, a truncation discussed 
extensively elsewhere [41, 47, 48, 50, 55-58]. This truncation requires to keep only terms 
up to one-loop, which include at least one ghost line. Therefore, also the gluon tree-level 
term can be dropped self-consistently. The equations then read 

D''t\p) = -ZsS^'V (2.5) 

+9r'' J -^,ip,{-q.P, q - p)D%{p - q)D'§ {q)rf^''sf {-p, q,p-q) 

Df-\p) = -gr'^l ^ip,D'J{q)D'G'{p + q)rf^'f''{-q,p + q,-p), 

where d is the dimensionality, Dq is the ghost propagator, D^i, is the gluon propagator, 
is the possibly finite ghost wave function renormalization, and p^'^"^'"^'^ jg f^n ghost- 
gluon vertex. The latter is undetermined at this level of the truncation, and will be set 
to its bare counterpart. For SU(2) and SU(3) this appears to be an even quantitatively 
good approximation [49, 59, 60], and is at least self-consistent for arbitrary gauge groups 
[56-58, 61]. If an additional color-tensor proportional to the symmetric color tensor d"'^'^ 
would appear, this would get lost due to the contraction with the antisymmetric tree-level 
vertex. Hence, at least the leading infrared part will have for neither propagator color off- 
diagonal elements. Only at two-loop order, and thus infrared and ultraviolet sub-leading, 
this could happen. The results on the lattice presented in section 5 and most accurate in 
this intermediate energy domain show that this is not the case. 
Taking the propagators to be of the form 

(p) = -5'^^^gP^2-2k ^2.6) 
DfAp) = -S^' - ^) Azp-'-'' (2.7) 



p. 

it is possible to absorb the remaining tree-level term in (2.5) by implementing the boundary 
condition [41] of an infrared divergent ghost dressing function, i. e., setting l/S = in the 
language of [20]. The integrals can then be performed analytically to yield 

= g'CAAlAzlG{^,t,d)p"(^-''^~'--'' 
= g^CAAlAzIz{^,d)p-^''-''^-''\ 

The expressions Iq and Iz are functions depending solely on the exponents k and t, and 
the structure of the underlying space-time manifold, symbolized by the dependence on d, 
and can be found, e. g., in [48]. Counting powers of momentum yields the relation 

4-d , , 

-2K = t + (2.8) 

The remaining consistency condition 

Iz{K,d) = lG{K,t{K),d) 
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implies that the exponents are only depending on the space-time manifold, but not on the 
gauge algebra: At this level, the fact that the propagators behave like a power-law and 
the value of the exponent is independent of the gauge algebra. A dependence on the gauge 
algebra can only be induced if the ghost-gluon vertex would be different from the tree-level 
version, and at least quantitatively dependent on the gauge algebra. However, as long as 
a solution to (2.8) exists with k > 1/2, this would only affect the numerical value of the 
exponent, but not the qualitative behavior of the propagators^. This yields the conjecture 
of the qualitative independence of the infrared behavior of the propagators of the gauge 
algebra. 

The pre-factors are, however, gauge-algebra dependent by virtue of the equation 



This gives the further prediction of a scaling relation with the adjoint Casimir in the far 
infrared 

p2+dj^aa2j^aa ^ (2.IO) 

In particular, in the 't Hooft limit [62] of su(A^) Yang-Mills theory the left-hand side of 
(2.10) must be a constant as a function of the number of colors N . Hence, Aq and Az must 
both scale in some way as a function of g^CA to ensure (2.10). However, if the 't Hooft- 
limit is taken, both Aq and Az must either compensate their respective scaling behavior, 
or must be individually independent of N. The latter behavior is indeed observed in the 
lattice calculations later, notably in two dimensions where a scaling behavior is manifest. 
This implies that the effective coupling, given by [39] 



will be essentially independent of the gauge algebra, and any dependence can only be 
introduced by the experimental input a(;U^), which cannot be fixed inside the theory. 

Note that the propagators are still non-trivial in the N ^ oo limit, and their properties 
are not obtainable by perturbative calculations in the planar limit. Nonetheless, the non- 
perturbative planar limit of the Dyson-Schwinger equations is sufficient to obtain this 
result. The failure of perturbation theory in the large A'^-limit is expected, as argued in 
the introduction. However, the fact that a non-perturbative planar limit of the DSEs 
is sufficient for the asymptotic infrared behavior is a non-trivial result. This actually 
generalizes to the DSEs for any Green's functions [55-58]. Note that also the all-order 
power-counting analysis [56, 57] of DSEs and the exact renormalization group equations 
still obey the gauge-algebra independence trivially. 

Hence, at this level neither k nor q'^CaAqAz receive any 1/A^ corrections in the large 
A'^-limit. Any dependence of these quantities on N therefore indicate A^-dependent vertex 
corrections. 

^Note that k. = 1/2 corresponds to an infrared finite rather than an infrared vanishing gluon propagator. 
This is not a quahtative different solution, when the characteristic of the solution is taken to be the relation 
(2.8), rather than the individual behavior of the propagators [56, 57]. 



g^CAAlAz 



1 1 




a(p2) = g^CAp'D^a^'D' 
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When moving to the photon-ghost case, the situation is less simple to access. In this 
case, at low momenta all diagrams contribute equally [45, 51, 53, 63, 64] to the effective 
gluon screening mass and the effective ghost wave-function dressing. In particular, the two- 
loop diagrams, which contain four-gluon vertices, potentially provide an equally important 
contribution. By this, further dependencies on further invariant tensors of the gauge algebra 
may appear. However, this is again not a qualitative effect, since the infrared behavior is 
not altered, but it can be quantitatively significant. However, the lattice results below 
suggest that such quantitative effects are comparatively small. 

3 Lattice formulation 

The standard Wilson action [65] 



is valid for any Lie-group, provided the link matrices U^, building the standard plaquette 
U^i, , are given in the fundamental representation of the group used to represent the gauge 
algebra. Np is then the dimension of the fundamental representation, and the sum is over 
all plaquettes. The bare coupling constant g is encoded in the constant (3 as 



where a is the lattice-spacing. In the implementation employed here, the fundamental 
representation for the SU(A^) groups has been taken from [66], and for G2 the Macfarlane 
representation [67] has been used. Np takes then the values 2 to 6 for SU(2) to SU(6), and 
7 for G2, respectively. 

A significant problem in comparing the results for various gauge groups are the poten- 
tially different scales [33]. This can be solved most easily by expressing all quantities by 
dimensionless ratios. An alternative is setting the string tension to the same value, here 
chosen conventionally to be (440 MeV)^. In this case all quantities are expressed effectively 
in units of the string tension. It is then possible to select also the same physical volume 
by choosing the same extension in lattice units and /3 such that for all gauge groups the 
physical volumes in units of the string tension agree. However, for d > 3, the asymptotic 
string tension in case of the gauge group G2 vanishes [32]. Nonetheless, an intermediate 
string tension, which is also used for setting the scale in case of the gauge group SU(A^)^, 
is non- vanishing [11, 69], and will be used here for this purpose. 

Unfortunately, even a direct measurement of the intermediate distance string tension 
is difficult and nontrivial in three and higher dimensions, because various corrections to 
the string tension are present. The values for SU(A^) have been taken from [70]^. For G2 
similar results in three dimensions are by now available in [37], and agree reasonably with 

^The asymptotic string-tension is the one in the Ai'-ality regime [68], which is not necessarily the one 
observed at short distances of a few fermi. 

^Interpolated and extrapolated in /3 and N, where necessary. 
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Table 1. The /3 values used in the simulations, together with the corresponding values of a and 
g for a string tension of (440 McV)^. N is the size of the largest available volume in lattice units. 
P is the value of the plaquette from this volume. Note that in two dimensions the same value of 
the string tension implies the same value of the plaquette. Of course, this would require precisely 
matching /3 values. Here, an agreement to two digits in a was taken to be sufficient. 



d 


Group 


P 


a [fml 

L J 


[GeV] 

L J 


g GeV— 


Ca 


g'^CA [GeV^-'^] 

tJ ^ L J 


N 


P 


2 


SU(2) 


10 


0.18 


1.1 


0.70 


2 


0.97 


150 


854185f2') 


2 


SU(2) 


38.7 


0.089 


2.2 


0.72 


2 


1.0 


150 


0.9614970(7) 


2 


SU(3) 


28 


0.18 


1.1 


0.51 


3 


0.78 


112 


0.85865(2) 


2 


SU(3) 


100 


0.089 


2.2 


0.54 


3 


0.89 


112 


0.96011(2) 


2 


SU(4) 


52.9 


0.18 


1.1 


0.43 


4 


0.73 


104 


0.85903(2) 


2 


SU(4) 


182 


0.089 


2.2 


0.46 


4 


0.86 


104 


0.95885(4) 


2 


SU(5) 


84.7 


0.18 


1.1 


0.38 


5 


0.71 


104 


0.85882(2) 


2 


SU(5) 


302 


0.089 


2.2 


0.40 


5 


0.81 


104 


0.960298(4) 


2 


SU(6) 


117 


0.18 


1.1 


0.35 


6 


0.74 


44 


0.85085(3) 


2 


SU(6) 


442 


0.089 


2.2 


0.36 


6 


0.79 


44 


0.960425(4) 


2 


G2 


50 


0.18 


1.1 


0.58 


2 


0.68 


56 


0.86010(4) 


2 


G2 


175 


0.089 


2.2 


0.63 


2 


0.78 


56 


0.959993(5) 


3 


SU(2) 


4.24 


0.17 


1.2 


1.0 


2 


2.1 


64 


744563fl) 


3 


SU(2) 


7.09 


0.094 


2.1 


1.1 


2 


2.3 


64 


0.8531134(7) 


3 


SU(3) 


10.7 


0.17 


1.2 


0.81 


3 


1.9 


56 


0.72480(2) 


3 


SU(3) 


18 


0.094 


2.1 


0.84 


3 


2.1 


56 


0.844316(9) 


3 


SU(4) 


19.85 


0.17 


1.2 


0.68 


4 


1.9 


40 


0.720095(7) 


3 


SU(4) 


34.2 


0.094 


2.1 


0.70 


4 


2.0 


40 


0.846161(4) 


3 


SU(5) 


31.45 


0.17 


1.2 


0.60 


5 


1.8 


32 


0.716233(8) 


3 


SU(5) 


54.3 


0.094 


2.1 


0.62 


5 


1.9 


32 


0.844750(5) 


3 


SU(6) 


45.7 


0.17 


1.2 


0.55 


6 


1.8 


20 


0.71469(3) 


3 


SU(6) 


79.3 


0.094 


2.1 


0.56 


6 


1.9 


20 


0.84491(1) 


3 


G2 


18.8 


0.17 


1.2 


0.93 


2 


1.7 


24 


0.72395(3) 


3 


G2 


32.9 


0.094 


2.1 


0.94 


2 


1.8 


24 


0.850734(8) 



the estimation procedure described in [33], which was used here. In two dimensions, the 
string tension is determined by direct measurements, which reproduces the known results 
for SU(2) [71], and supports the method of [33] for G2. 

Still, a quantitative comparison between the different gauge groups makes only sense 
for inherently dimensionless quantities, like the infrared exponents, as the value of the 
string tension is arbitrarily set to the same value. Hence, only qualitative statements for 
dimensionful quantities are sensible. Note that this might even affect the determination of 
the volumes. The selected values for /3, and the corresponding values of a and g, are given 
in table 1. Since the coupling is the only unique dimensionful quantity in these theories, it 
would be equally well valid to set it to some fixed value, say 1 GeV, and determine by this 
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Table 2. The number of configurations for the various systems investigated. The number of 
configurations is given for the various gauge groups. Swe. is the number of sweeps between two 
consecutive measurements, and Ther. is the number of initial thermahzation sweeps. Note that 
always multiple independent runs have been performed to obtain the final statistics. The latter was 
selected to obtain the ghost exponent with a statistical accuracy of 10% at the Icr- level, compu- 
tational resources permitting. High and low in the /3 column corresponds to the higher and lower 
P value for the different gauge groups. An empty field indicates that no measurements have been 
performed for this gauge group at this volume and (3 value. V is the physical volume and the 
lattice extension. The data for SU(2) are taken from [21]. 



d 


P 


N 


1 

Vd [fm] 


Ther. 


Swe. 


SU(2) 


SU(3) 


SU(4) 


SU(5) 


SU(6) 


G2 


2 


High 


4 


0.36 


140 


14 


1045 


935 


777 


620 


465 


620 


2 


Low 


4 


0.72 


140 


14 


1045 


935 


876 


620 


459 


620 


2 


High 


10 


0.89 


200 


20 


1589 


1087 


840 


567 


483 


543 


2 


High 


16 


1.4 


260 


26 


1017 


895 


735 


641 


451 


597 


2 


Low 


10 


1.8 


200 


20 


2365 


1169 


848 


567 


519 


584 


2 


High 


20 


1.8 


300 


30 


1094 


895 


744 


642 


421 


562 


2 


High 


26 


2.3 


360 


36 


1015 


903 


807 


629 


442 


519 


2 


Low 


16 


2.9 


260 


26 


1017 


898 


734 


635 


403 


610 


2 


High 


34 


3.0 


440 


44 


1038 


879 


712 


620 


496 


566 


2 


Low 


20 


3.6 


300 


30 


1041 


896 


788 


626 


433 


570 


2 


High 


44 


3.9 


540 


54 


1060 


931 


769 


581 


147 


265 


2 


Low 


26 


4.7 


360 


36 


1009 


903 


774 


581 




554 


2 


High 


56 


5.0 


660 


66 


1054 


912 


718 


359 




68 


2 


Low 


34 


6.1 


440 


44 


1038 


1756 


795 


575 


380 


552 


2 


High 


68 


6.1 


780 


78 


1011 


918 


516 


25 






2 


High 


80 


7.1 


900 


90 


1113 


358 


94 


51 






2 


Low 


44 


7.9 


540 


54 


1059 


9954 


1476 


751 


59 


318 


2 


High 


104 


9.3 


1140 


114 


1045 


30 


41 


25 






2 


Low 


56 


10 


660 


66 


6564 


4611 


853 


491 




19 


2 


High 


112 


10 


1220 


122 


3524 


10 










2 


Low 


68 


12 


780 


78 


1011 


1694 


342 


87 






2 


High 


136 


12 


1460 


146 


2179 












2 


Low 


80 


14 


900 


90 


1048 


274 


115 


38 






2 


High 


150 


13 


1600 


160 


1066 












2 


Low 


104 


19 


1140 


114 


1032 


50 


48 


28 






2 


Low 


112 


20 


1220 


122 


2250 


45 










2 


Low 


136 


24 


1460 


146 


1981 












2 


Low 


150 


27 


1600 


160 


2238 













condition a. As the values in table 1 indicate, this would modify the scales by up to 50%. 

Configurations are then obtained by a mix of heat-bath and overrelaxation sweeps. 
This has been done for SU(2) as described in [72]. For SU(A^ > 3), heat-bath updates 
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Table 3. As table 2, but for three dimensions. The SU(2) data is again taken from [21]. 



d 




N 


[fm] 


Ther. 


Swe. 


SU(2) 


SU(3) 


SU(4) 


SU(5) 


SU(6) 


G2 


3 


High 


4 


0.36 


240 


24 


1070 


982 


832 


574 


465 


574 


3 


Low 


4 


0.72 


240 


24 


1070 


975 


771 


574 


405 


574 


3 


High 


8 


0.72 


280 


28 


1133 


967 


707 


684 


230 


520 


3 


High 


12 


1.1 


320 


32 


1128 


932 


730 


624 


410 


531 


3 


Low 


8 


1.4 


280 


28 


1070 


967 


701 


684 


460 


596 


3 


High 


16 


1.4 


360 


36 


1144 


982 


774 


572 


435 


648 


3 


High 


20 


1.8 


400 


40 


1100 


931 


706 


570 


46 


285 


3 


Low 


12 


2.2 


320 


32 


1128 


930 


724 


570 


422 


507 


3 


High 


24 


2.2 


440 


44 


1084 


939 


970 


307 




68 


3 


Low 


16 


2.9 


360 


36 


1048 


985 


786 


572 


229 


616 


3 


High 


32 


2.9 


520 


52 


960 


935 


270 


60 






3 


Low 


20 


3.6 


400 


40 


1217 


979 


734 


617 


33 


172 


3 


High 


40 


3.6 


600 


60 


1065 


102 


74 








3 


Low 


24 


4.3 


440 


44 


1238 


867 


719 


276 




35 


3 


High 


48 


4.3 


680 


68 


1095 


100 










3 


High 


56 


5.0 


760 


76 


1078 


16 










3 


Low 


32 


5.8 


520 


52 


899 


665 


193 


62 






3 


High 


64 


5.8 


840 


84 


1062 












3 


Low 


40 


7.2 


600 


60 


1040 


156 


75 








3 


Low 


48 


8.6 


680 


68 


1047 


74 










3 


Low 


56 


10 


760 


76 


1053 


17 










3 


Low 


64 


12 


840 


84 


1033 













have been performed using the Cabibbo-Marinari method [73] with 3, 6, 10, and 15 SU(2) 
subgroups for = 3, 4, 5, and 6, respectively, according to [33]. Between two heat- 
bath sweeps five overrelaxation sweeps have been performed. In the case of G2, heat-bath 
updates have been performed according to the method presented in [10] and detailed in 
[11], and overrelaxation sweeps according to [33]. This has been done in both, two and 
three dimensions. The number of configurations are given in table 2 and 3 for two and 
three dimensions, respectively. 

Once an equilibrated configuration is obtained, it is still necessary to fix it to the 
Landau gauge. This is done by minimizing the functional (see, e.g., [19]) 

£ = -^mTUf,, (3.1) 

X,fl 

where {Ufj_{x)} is a thermalized lattice configuration, and the sum is over all links. This 
leads to the so-called minimal [19], or average- [20], Landau gauge. In principle, for 
comparison to the results of section 2.3, it is necessary to fix to the non-perturbative 
version of the Landau gauge which exhibits scaling, if such a gauge indeed exists [42]. For 
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the volumes available here for SU(A^ > 2), this appears to be just a quantitative effect in 
the SU(2) case for the actual propagators [20, 21], and only slightly modifies the volume 
dependence. For the purpose here, investigating the changes with the gauge group, these 
effects are sub-leading at the present volumes, and therefore the much cheaper minimal 
Landau gauge is sufficient. However, it is not known yet, whether the differences between 
the various Landau gauges is stronger for gauge groups different from SU(2). Studies in 
four dimensions with SU(3) do not suggest so [74]. Hence, this will be assumed henceforth. 
In particular, it will be assumed that a qualitatively similar result for the gauge-group 
dependence in minimal Landau gauge for the volumes investigated here implies the same 
qualitative behavior in other Landau gauges as well. 

It is then sufficient to fix to the minimal Landau gauge. For SU(2), this is done using 
stochastic overrelaxation with adaptive parameter adjustment [72]. For the other gauge 
groups, only the standard overrelaxation step has to be modified compared to the one in 
SU(2). This can be done by overrelaxing all SU(2) subgroups, in the same manner as for 
overrelaxation sweeps during the generation of configurations, in particular the same num- 
ber of subgroups. For SU(3) the procedure is described in [75], and can straightforwardly 
be generalized to higher values of N. For G2 this is described in [33]. The quality of the 
gauge-fixing is monitored using the quantity eg, defined as [19] 



where {g{x)} represents the gauge transformation applied on the link variables U^{x), the 
symbol + indicates Hermitian conjugation, is the lattice side of the symmetric hypercube, 
d is the space-time dimensionality, is a positive unit vector in the fi direction and Tc are 
the generators of the algebra. This quantity is a more reliable measure of the gauge-fixing 
quality than just the transversality itself [19]. Furthermore, it is found that the same limit 
of eg corresponds to a much better fulfillment of the transversality condition with increasing 
number of generators. Therefore, the restriction on eg for achieving the gauge-fixing can 
be taken somewhat lower for G2 and SU(A^ > 2) than for SU(2), corresponding still to a 
better level of transversality on the average. The limits adopted here are 10"^^ for SU(2) 
and 10~^^ otherwise. 

4 Propagators 

The determination of the propagators proceeds in the same way for all gauge groups. 
The gluon propagator is given by the correlation function 




(3.2) 




[9{x)Uf,{x)g{x + 6/,)+ - g{x + ei,)Ui,{x)'^ g{x)+] 




(4.1) 
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with the momentum-space lattice gluon field defined as 

Here the components p^j of p have the integer values —N/2 + 1 , . . . , N/2 . 

After contracting (4.1) with a transverse projector and a unit matrix in color space, 
the scalar part of the gluon propagator is given by 

/J, a 

where 5ftAJ^(p) and 9^4^ (p) are, respectively, the real and the imaginary part of ^JJ(p) and 
the normalization M is given by dNg for p > and by {d — l)Ng for p = 0. Ng is the 
number of generators and thus the number of gluons, for a given gauge group. Hence, 
Ng = iV^ — 1 for SU(A^) and Ng = 14 for G2. The gluon propagator is thereby by definition 
positive semi-definite. 

The components of the physical momenta are given by 

P^ = 2sin^. (4.4) 

Results will be presented as a function of the magnitude of the physical momentum 
p = \P\/a (in GeV). Note that the continuum gluon propagator is obtained by the product 
/3 a?D{k), since the lattice quantity ^/p/a/^ A'^^{x) yields the continuum quantity ^^(x) in 
the formal continuum limit a — > 0. In the same limit, (3 a'^^'^A'^{p) converges to the con- 
tinuum momentum-space gluon field j4^(p). Thus, for any dimension d, the lattice quantity 
j3 a^Dik) converges to the continuum gluon propagator in momentum space, independently 
of the gauge group. 

The ghost propagator is given by 

D^ciP) = ^ < {M-^T\p) >, (4.5) 

where M°'^{x,y) is the Faddeev-Popov operator, defined in the continuum as 

- d^Df = 6{x - y){-dH-' + gr'^d^A^) . (4.6) 

On the lattice in Landau gauge this operator is a matrix defined by its action on a scalar 
function w^(x) as [76] 



M(y, xT'uobix) = c\Y, {G''\x)uj,{x) + {x)ujbix + e^) + <(x)a;b(x - e, 

\ X II 



)) 



C^'ix) = Y^r{{r^^r'}{U,{x) + U,{x-e,))) 

Afix) = -2tTiT^T'U,ix)) 
Bf{x) = -2tr(TV^+(x-e^)), 
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where c is a constant depending on the normahzation of the generators of the gauge algebra. 
The evaluation of the Fourier transform of the inverse operator 

with p = —q requires a matrix inversion, which has been performed for all gauge groups 
using the point source d°''^{5xo — [77]. Independent of the gauge group the Faddeev- 

Popov operator M°'^[x,y) is symmetric and positive. Hence, in practice the matrix inver- 
sion has been performed using a conjugate-gradient method [72]. Fortunately, the larger 
the number of generators the smaller is the extensive statistical noise induced by using a 
point source instead of a plain wave source. However, better accuracy is required in the 
inversion process such that the spectrum of the Faddeev-Popov operator is positive with 
an increasing number of generators. Otherwise negative eigenvalues are encountered as nu- 
merical artifacts, which become positive when increasing the accuracy [78]. This procedure 
is ambiguous with respect to the sign of the resulting propagator (or of the eigenvalues) 
[72]. Thus, the sign has to be assigned by hand, and is fixed by comparing the propagator 
at large momenta to perturbation theory. 

Finally, the color- averaged ghost propagator is then defined as 

Dcip) = -^D^G^. (4.8) 

One particular question is whether it is sufficient to investigate the color-averaged propa- 
gators (4.3) and (4.8). For SU(2), this is the case [49, 72], and is not surprising due to the 
vanishing of the symmetric structure constants. For other gauge groups, there is no reason 
why this should be the case. In the next section, it will be shown that this justified. 

5 Results 

5.1 Propagators 

The initial check is whether the color-off-diagonal propagator elements indeed vanish, as 
this was a necessary precondition for the analysis of section 2.3, and is a general assumption 
in functional studies. The ffi'st results are therefore for the color off-diagonal, but Lorentz- 
structure averaged, gluon propagator in figure 1 in two dimensions and in figure 2 for three 
dimensions, for the various gauge groups. In both two and three dimensions and for all 
gauge groups the results are, within statistical errors, in agreement with a zero result. The 
gluon propagator is thus color-diagonal, at the very least for the lattice settings studied 
here. 

The corresponding result for the ghost propagator is shown in figure 3 for two dimen- 
sions and in figure 4 for three dimensions. Once more, for all gauge groups the off-diagonal 
elements are consistent with zero, in agreement with the arguments of section 2.3. 

This supports corresponding assumptions in functional calculations of all types. Note 
that this is not necessarily implying a simple color-structure on the level of the vertices 
[79, 80]. To some extent, this is an expected result, as off-diagonal color elements could 
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Figure 1. The ofF-diagonal gluon propagator in two dimensions for tlie different gauge groups. 
Always the components D^^ (full symbols) and Z?^^? (open symbols) are shown. 



be indicative of the breaking of the residual global color symmetry, which would be unex- 
pected. However, note that it is currently unclear whether in the minimal Landau gauge 
the corresponding charge can be well-defined for a photon-ghost behavior [81]. 
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Figure 2. The ofF-diagonal gluon propagator in tiiree dimensions for tlie different gauge groups. 
Always the components D^^ (full symbols) and Z?^^? (open symbols) are shown. 



The results for the color-diagonal gluon propagator is shown in figure 5 for two dimen- 
sions and in figure 6 for three dimensions. 

In two dimensions, for all gauge groups, the gluon propagator is infrared suppressed. 
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Figure 3. Same as in figure 1, but for tlie gliost propagator. 



compared to tree-level. In fact, for sufficiently large volumes it is also stronger suppressed 
than the one of a massive particle, as expected from section 2.3 and the previous results 
for SU(2) alone [25, 49]. Furthermore, the relevant scales, expressed in units of the string 
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Figure 4. Same as in figure 2, but for tlie gliost propagator. 



tension, turn out to be essentially the same for all gauge groups: The maximum occurs at 
about half a GeV, and the height of the maximum is between two and three inverse GeV^. 
Also, in all cases the effects of the violation of rotational symmetry are rather small. The 
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Figure 5. The gluon propagator in two dimensions for the various gauge groups. Error bars are 
partly smaller than the symbol size. Various momentum configurations are shown, see [72] for 
details. 



infrared suppression observed is in all cases also sufficient to cure tlie infrared problems 
encountered in perturbation theory, making all integrals in the Dyson- Schwinger equations, 
or elsewhere, well-defined. 

In three dimensions, only in the SU(6) (and possibly the G2) case an explicit maximum 
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Figure 6. Same as in figure 5, but for three dimensions. 



is not visible, although the finite- volume effects for all gauge groups demonstrate explicitly 
that D(0) has still not settled on its infinite- volume value. Hence, the existence of a 
maximum at larger volumes for the group SU(6) appears likely. Furthermore, for the same 
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Figure 7. The ghost dressing function for the various gauge groups in two dimensions. Error bars 
are partly smaller than the symbol size. Various momentum configurations are shown, see [72] for 
details. 



volume for all other gauge groups also no maximum is visible. Hence, the qualitative 
turn-over structure seems to agree for all gauge groups in three dimensions. 
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Figure 8. Same as in figure 7, but for three dimensions. 



The results for the ghost dressing function are shown in figure 7 for two dimensions 
and in figure 8 for three dimensions. 

In two dimensions again there is no quahtative, and httle quantitative, difference be- 
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Table 4. Fit results for the ultraviolet fits. Errors are statistical only. 



Propagator 


Group 


d 


a 


a 

Ca 


a 

Caq^ 


d 


a 


a 

Ca 


a 

CAd^ 


a 
OLD 


Gluon 


SU(2) 


2 


-0.43(11) 


-0.22 


-0.43 


3 


-0.59(3) 


-0.30 


-0.26 


1.5 


Gluon 


SU(3) 


2 


-0.46(12) 


-0.15 


-0.52 


3 


-0.58(6) 


-0.19 


-0.28 


1.6 


Gluon 


SU(4) 


2 


-0.47(6) 


-0.12 


-0.55 


3 


-0.56(5) 
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-0.28 


1.6 


Gluon 


SU(5) 


2 


-0.43(6) 
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-0.53 


3 


-0.56(4) 
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1.7 
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3 
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-0.28 


1.7 
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3 
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tween all gauge groups. The ghost dressing function is in all cases infrared enhanced, and 
compatible with being infrared divergent. Also, the effects of violation of rotational in- 
variance are once more rather small. Furthermore, the ghost dressing function is for all 
gauge groups monotonous for the complete momentum regime. Thus, as in case of the 
gluon propagator, the similarity between the various gauge groups holds also beyond the 
asymptotic regimes discussed in section 2. 

In three dimensions, in all cases, the dressing function is infrared enhanced. In case 
of SU(2), substantial evidence exists that this enhancement is not a divergence, and the 
ghost dressing function becomes finite at very small momenta [24], despite its apparent 
strength. Since the propagators for the other gauge groups are similar to the one of SU(2) 
at equivalent volumes, though less strongly enhanced, they will likely become also finite at 
larger volumes. Nonetheless, the main result here is that these ghost dressing function are 
qualitatively very similar, and the influence of the gauge algebra on the ghost propagator 
at these volumes is therefore only quantitative. 

5.2 Asymptotic behavior 
5.2.1 Ultraviolet behavior 

The first interesting question to be established is whether the universality at large momenta 
obtained perturbatively in higher dimensions also holds in two dimensions. For this it is 
necessary to obtain the leading perturbative corrections to the tree-level value. This is 
done by fitting the data for both propagators above 2 GeV and for the volumes used in the 
figures of section 5 and the larger /? value for momenta along the diagonal with the form 

p2i^(p) = (5.1) 

The results of these fits are given in table 4. The results for the gluon propagator are found 
to cluster around a value of c = ajig^CA) of about 1/2. The values are always negative. 
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Figure 9. The effective gluon (left) and ghost (right) exponents in two dimensions. Full circles are 
SU(2), full squares SU(3), full triangles SU(4), full upside-down triangles SU(5), open circles SU(6) 
and open squares G2. 



demonstrating the existence of a Landau pole also in two dimensions. Using instead of 2 
GeV 4 GeV as the starting point for the fit results in a slight increase in the value of a, 
but at the expense of a larger statistical error. A similar results is found for the ghost 
propagator, though the value is now closer to 3/10. Altogether, the results show that 
the perturbative behavior in two dimensions follows the expected pattern of gauge algebra 
dependence. 

To compare how well this approximates the leading-order behavior also the results 
for three dimensions are given and compared to the leading-order result. For the gluon 
propagator, the deviation is sizeable^, but still not dominant. The agreement is slightly 
but not drastically improved when starting to fit at 4 GeV instead. 

For the ghost propagator, the leading-order approximation essentially breaks down at 
2 GeV, and is not much better at 4 GeV. The results, which therefore include higher-order 
corrections, still show an approximate scaling proportional to q^Ca with the gauge-group, 
showing that even sub-leading corrections are pre-dominantly of this type, rather than 
having more complex dependencies on the gauge algebra. 

Note that despite these differences the propagators above 2 GeV are almost completely 
dominated by their tree-level behavior. Even with the enhancement in the ghost case, the 
leading-order corrections at 2 GeV amount only to a 15-20% effect. Without it, it is only 
a 5% effect. 

5.2.2 Infrared behavior in tw^o dimensions 

As shown in figures 5 and 7, the propagators exhibit for all gauge groups a behavior 
which is similar to the SU(2) case in approximately similar volumes. Therefore, a fit with 

*Some of the deviations could possibly be accommodated by a tadpole correction [82], which amount to 
about a ten-precent effect at these values of jS. 
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the ansatze (2.6) and (2.7) can be performed, as described in [49]. This yields volume- 
dependent effective exponents. They are shown in figure 9 for the various gauge groups. 

The curves for the gluon exponent are rather similar, and the small differences could 
easily be dominated by the uncertainties in the comparison of the scales. In particular, the 
value of the exponent at large volumes appears to be rather insensitive to the gauge group. 

The same observations pertains to the ghost exponent. In general, the value of the 
exponent is in agreement between all gauge groups, and seems to be consistent with the 
same value in the infinite- volume limit. With the rather larger statistical errors, this cannot 
be finally settled 

Note that in both cases these evolution statements can only be made under the as- 
sumption that the assigned volume scale is indeed in rough agreement between all gauge 
groups. Therefore, the most important statement is that the infinite-volume limit seems 
to be in agreement, as expected from the discussion of section 2.3. This also implies that 
gauge-group-dependent corrections to the low-momentum behavior of the ghost-gluon ver- 
tex are likely small and of quantitative nature only. This is in agreement with investigations 
of the ghost-gluon vertex for SU(2) and SU(3) in four dimensions using lattice gauge theory 
[59, 60], as well as with functional studies [61]. 

Although not in all cases the fitted gluon exponents have yet reached a stable value as 
for SU(2), a maximum is already visible in all cases. If the gluon propagator is to vanish 
in the infinite- volume limit in two dimensions, then the value of it at zero momentum 
is predicted to behave as a power of inverse volume [83]. This is indeed the case for all 
volumes investigated here, as can be seen in figure 10. In fact, for all gauge groups the 
gluon propagator vanishes like the same inverse power of V. This inverse power is of the 
same size as it would be expected on the basis of finite volume corrections obtained in 
functional calculations [49, 83]. 
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Figure 11. The coupling (2.11) for tlie various gauge groups in two dimensions. Always the results 
for the largest volume with the low /3 value (crosses) and the high /3 value (full circles) are shown. 
In addition, to assess discretization effects, also the result for the lower /3- value, which has a volume 
most closest to the result for the higher /3 value shown, is presented as well. 
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The final quantity to be evaluated is the effective running coupling, which is given 
in (2.11). The results for the various gauge groups are shown in figure 11. Not for all 
gauge groups sufficiently large volumes are available to identify the existence of an infrared 
plateau, as is possible for SU(2). However, the approximate value is in all cases about the 
same, and thus, if at all, there is only a weak dependence of the effective coupling in the 
infrared on the gauge group. 

5.2.3 Infrared behavior in three dimensions 

In three dimensions in minimal Landau gauge the behavior of the propagators change to a 
photon-ghost behavior at a momentum of order a few dozen MeV [24, 25], though still at 
momenta much smaller than in four dimensions [23-25]. Ultimately, therefore, in minimal 
Landau gauge the gluon propagator resembles the one of a screened particle and the ghost 
propagator becomes photon-like. Still, this critical momentum is lower than accessible 
in most cases here. Only for SU(2), this momentum is just so reached. But since such 
momenta are currently almost not accessible for most gauge groups, and the intermediate 
low-momentum behavior appears to be power-like for all groups, this will be investigated 
here. Since in this intermediate momentum range a scaling-like behavior is expected [41], 
the analysis of section 2.3 can be applied with an infrared cutoff [83]. As a consequence, 
the behavior of the propagators in this intermediate range should again be universal, i. e., 
independent of the gauge algebra up to trivial 't Hooft scaling. This will be checked here. 

Taking thus the low-momentum forms (2.6) and (2.7) as fit ansatze, the effective, 
volume-dependent exponents can be determined. These are plotted in figure 12. At least 
for SU(2) and SU(3) the gluon exponent becomes smaller than -1 for the given volume, i. 
e., at least for some momentum range the gluon propagator decreases faster than would 
be expected from a screened particle behavior. Of course, in much larger volumes this 
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Figure 13. The gluon propagator at zero momentum in tliree dimensions. Full circles arc SU(2), 
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open squares G2. 



exponent increases again towards -1. Furthermore the volume-dependence of the effective 
exponents is once more the same for all gauge groups, as is expected from the presence of 
the scaling window. However, quantitatively, there is some difference in the ghost case for 
the gauge groups. In addition, the gluon propagator at zero momentum should decrease 
with volume for a certain range of volumes, like a power-law. This is indeed the case, as 
can be seen in figure 13. 

Finally, it is worthwhile to also investigate the effective coupling. The results are 
shown in figure 14. Again, the behavior is also quantitatively rather similar. In all cases, 
the coupling is infrared suppressed, as is expected in any case. However, the maximum can 
be compared to (2.9). Since the highest value reached is more or less similar this implies 
that the maximum scales approximately like (2.9). The difference between SU(2) and the 
remaining groups maybe connected with the fact that the peak is here at about 400 instead 
of 500-600 MeV. There, again, some scale uncertainty may be involved. 

Hence, also in three dimensions, the asymptotic behavior at low momentum follows in 
an intermediate volume window a scaling-type behavior, as is expected independently of 
the far infrared behavior. Even without investing this information, the result is that the 
same pattern emerges for all gauge groups, and thus even quantitatively the effects of the 
different gauge algebras are small. Furthermore, the behavior is in all cases not compatible 
with perturbation theory at low momenta, nor do non-perturbative contributions seem to 
diminish with A^. Thus it is likely that even in the limit of — > 00 non-perturbative effects 
dominate the infrared behavior of the correlation functions, at least for momenta below 

Of course, this is exactly as it would be expected when topological configurations, 
contributing with an essential singularity in the coupling constant, dominate the low- 
momentum behavior. There is evidence for this [84-88], so this behavior is not in disagree- 
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Figure 14. Same as figure 11, but for three dimensions, 
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6 Summary 



Summarizing, the propagator in two and three dimensions are quahtatively, and actuahy 
also quantitatively when measured in units of the string tension, very similar, for the gauge 
groups SU(2), SU(3), SU(4), SU(5), SU(6) and G2, and thus for the corresponding gauge 
algebras. In particular, no pronounced dependence on N is observed for the approximately 
fixed value of q^Ca used here. Hence, the dominating contribution in the mid-momentum 
regime and at low momenta are of order 0(1) in terms of N counting. In particular, for 
all N the propagators show a behavior which is distinctively non-perturbative. 

Furthermore, the results are compatible with predictions for the gauge-group depen- 
dence expected from functional calculations. In particular, the deep infrared remains the 
same for any gauge group. 

These results emphasize that for properties of Yang-Mills theories on the level of gluonic 
correlation functions the specific gauge algebra is qualitatively, and to some extent even 
quantitatively, rather irrelevant. Only when matter is coupled to the Yang- Mills field 
dynamically the gauge algebra structure becomes quite relevant (see e. g. [89]). Of course, 
even for pure Yang- Mills theory the gauge algebra is important for some quantities, e. g., 
the order of the finite-temperature phase transition. 
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